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We report an experimental investigation of the longitudinal space-time cross-correlation 
function of the velocity field, C(r, r), in a cylindrical turbulent Rayleigh-Benard convec- 
tion cell using the particle image velocimetry (PIV) technique. We show that while the 
Taylor's frozen-flow hypothesis does not hold in turbulent thermal convection, the re- 



cent elliptic model advanced for turbulent shear flows (|He fc Zh ang 2006) is valid for the 
present velocity field for all over the cell, i.e., the isocorrelation contours of the measured 
C(r, t) have a shape of elliptical curves and hence C(r, r) can be related to C(te, 0) via 
r% = (r — Ut) 2 + V 2 t 2 with U and V being two characteristic velocities. We further 
show that the fitted U is proportional to the mean velocity of the flow, but the values of 
V are larger than the theoretical predictions. Specifically, we focus on two representative 
regions in the cell: the region near the cell sidewall and the cell's central region. It is 
found that U and V are approximately the same near the sidewall, while U ~ at cell 
center. 

Key Words: Rayleigh-Benard Convection, velocity space-time correlations, Taylor's 
frozen-turbulence hypothesis 



1. Introduction 

Turbulent flows contain eddies with various scales. Turbulent kinetic energy is trans- 
ferred from eddies with the largest scale of turbulence, L, at which energy is injected into 
the turbulence system, to eddies with the smallest scale of turbulence, ~q, at which energy 
is dissipated by fluid viscosity. Such cascade processes are usually characterized by the 
velocity structure functions, S p (r) — {\5 r v\ p ), defined as moments of velocity increments 
ov er a space separatio n r, where (• • •) denotes a time average. Since the pioneering work 
of iKolmogorovl (jl94l[) , various theories and models have been put forwards to predict 
the scaling behaviors of the veloc ity structure functions in the so-called inertia ! range 
rj << r « L (see, for reviews, iFrisch. IJ.I Hiffli IS reenivasan fc Antonial 119971 ) . From 



experimental aspects, velocity measurements are usually carried out at a single fixed lo- 
cation, based on which time series of fluctuating velocities are obtained and the velocity 
structure functions, S p (t) = (\S T v\ p ), are calculated as moments of velocity increments 
over a time separation r. To relate the properties of the experimentally measured S p (t) in 
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time domain to theoretical predi ctions of S v ( r) in space domain, one needs to invoke the 
Taylor's frozen- flow hypothesis ( Taylor 19381 ). The validity of such hypothesis demands 



low turbulent intensity and weak shear rates ( Lumlevl 1965). Howeve r, these conditions 



are not always met by actual flows of interest (jPinton fc Labbdll994l ). 

Another quantity that can also be used to characterize the cascade processes is the 
velocity space-time correlation function, defined as 

C(r,T)= {v{x+ . r : t + T ) v{x ^ , (1.1) 

where v(x) is one of the components of the velocity vector at position x and v rms (x) 
is the root-mean-square (r.m.s.) velocity at x. For simplicity, we consider here only the 
situation that r is in the direction of v(x), i.e. the longitudinal velocity correlations. 
[Studies of other kinds of velocity correlations, such as those between t he wall-normal 
veloci ty component and the streamwise component, could be found in Ijachens et al 



(2006).] Again, the velocity temporal auto-correlation function C(0, r), based on the 



pointwise measurements, is the quantity most often studied in experiments and Taylor's 
frozen-flow hypothesis is usually used to translate C(0,r) in time domain to C(r, 0) in 
space domain. When Taylor's hypothesis holds, we have C(r, r) = C(rT,0) with 

r T = r- U t, (1.2) 

where Uq is the mean velocity of the flow. Such a relation implies that C (r, r) would 
keep constant at increasing r and r once the value of r — Uqt remains constant, which 
violates the basic properties of the correlations that C(r, r) decays to zero at sufficient 
large separations. Therefore, Taylor hypothesis holds only for the limited ranges of r and 

T. 

Recently, based on a second order approximation to C(r, r), Irle fc Zhangl ( 20061 ) ad- 



vanced an elliptic model for turbulent shear flows and proposed that C(r, r) could be 
related to C(te,0) via 

r% = {r-Urf + V 2 t 2 , (1.3) 

where U is a characteristic convection velocity proportional to the mean velocity Uq and 
V is a c haracteristic veloc ity associated with the r.m.s. velocity and the shear-induced 
velocity ( Zhao fc Hel 12009) . Specifically, when V vanishes, (|1.3p is degener ated to the 



Taylo r's hypothesis (II. 2p . while Kraichnan's sweeping- velocity hypothesis ( Kraichnanl 



Il964l ) is obtained if U vanishes. The He's elliptic model provides a useful tool for a large 
set of flow systems where the Taylor's hypothesis does not hold and hence the validity 
of the model in various practical flows needs to be tested. 

In this paper, we want to validate the elliptic model in turbulent convection, an impor- 
tant class of turbulent flows that play an central role in many natural and engineering 
processes. The flow at hand is turbulent Rayleigh-Benard (RB) convection, i.e. the con- 
vective motion of an enclosed fluid layer heated from below and cooled from above, which 
has received tremendous attention during the past few decades (Ahlers, Grossmann & 
Lohse 2009; Lohse & Xia 2010). Although it has long been recognize d that the condition s 
for Taylor's hypothesis are often not met in the system (see, e.g., IShang fc Xial 2001 ). 



single-point or time-domain measurements are employed by most experimental investiga- 
tors for the studies of turbulent cascade processes of both the velocity and temperature 
fields (see the recent review paper, Lohse & Xia 2010, and references therein). To trans- 
late correctly the quantities measured in time domain to those in space domain, it is thus 
essential to validate the elliptic model in turbulent RB convection. 

Recently, He, He & Tong (2010) verified indirectly the elliptic model via the local 
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temperature data as in the bulk region of turbulent RB convection temperature behaves 
as a passive scalar (Calzavarini, Toschi & Tripiccione 2002; Xi et al. 2009), which is driven 
by the velocity field via a linear equation. The authors showed that the elliptic relation 
(|1.3[) is valid for the temperature space-time correlations measured in the cell's sidewall 
region, but at the cell center they did not observe the predicted relation (jl.3|) for the 
temperature field. However, a passive ad ditive may disp lay characteristics so different 
from those of the advecting velocity field ( Warhaftll2000i ). In addition, the temperature 
field is nearly homogeneous in the cell's central region. The lack of temperature contrast 
would make the local temperature not follow the behaviors of the veloc ity field, e . g. the 
local velocity flu ctuations show strong oscillations at the cell center (|Qiu et al\ 12004 ; 



Zhou et all 2009), whi le the oscillation is absent for temperature measured at the same 



location ( Qiu fc Tone! 2002 ). and hence it is not surprising that (|1.3[) does not hold 
for temperature in the cell's central region. Therefore, it is highly desirable to verify the 
elliptic model directly via the velocity field, which is the object of the present experimental 
investigation. 

The remainder of this paper is organized as follows. We give detailed descriptions of the 
experimental apparatus and conditions and the measuring technique in §2. Experimental 
results are presented and analyzed in §3, which is divided into three parts. In §§3.1 and 
3.2, we study the properties of longitudinal space-time correlation functions of the vertical 
velocity near the cell sidewall and at the cell center, respectively. Section 3.3 presents 
results of longitudinal space-time correlations of the vertical velocity along the cell's 
diameter at the middle height of the cell and of longitudinal space-time correlations 
of the horizontal velocity along the cell's central vertical axis. The properties of the 
characteristic velocities U and V are also investigated in detail in §3.3. We summarize 
our findings and conclude in §4. 



2. Experimental setup and parameters 

The convection cell is similar to that us ed in previous experiments (|Sun et al.ll2005al ) 



but has a different size (IZhou fc Xiall2010h . It is a vertical cylinder of height H — 50 cm 
and inner diameter D = 50 cm and hence of unit aspect ratio. Dcionized and degassed 
water was used as the convecting fluid. The cell's sidewall is a plexiglas tube of 5 mm 
in wall thickness and a square-shaped jacket made of flat plexiglas plates and filled with 
water is fitted round the sidewall, which greatly reduced the distortion effect to the PIV 
images caused by the curvature of the cylindrical sidewall. The top and bottom plates 
are made of pure copper with nickel-plated fluid-contact surfaces. The thickness of the 
top plate is 3 cm and that of the bottom one is 1.5 cm. Four spiral channels of 1.2 cm in 
width and 1.5 cm in depth are machined into the top plate and the separation between 
adjacent channels is 1.1 cm. The channels start from the center and end near the edge of 
the plate. A silicon rubber sheet and a Plexiglas plate are fixed on the top to form the 
cover and also to prevent interflow between the channels. Each channel is connected to a 
separate refrigerated circulator (Polyscience 9712) that has a temperature stability of 0.01 
°G. The channels and the circulators are connected such that the incoming cooler fluid 
and the outgoing warmer fluid in adjacent channels always flow in opposite directions. 
Four quarter-circular Kapton film heaters, connected in parallel to a dc power supply 
(Xantrex XDC 300-20) with 99.99% long-term stability, are sandwiched to the back side 
of the bottom plate to provide constant and uniform heating. Therefore, the experiments 
were conducted under constant heating of the bottom plate while maintaining a constant 
temperature at the top plate. Eight thermistors are embedded beneath the fluid-contact 
surface of each conducting plate, equally spaced azimuthally at about one-third radius 
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Figure 1. Sketch of the convection cell, the Cartesian coordinates used in the experiment, and 
the time-averaged vector map of the whole velocity field measured in the vertical plane of the 
large-scale circulation at Ra — 9.5 x 10 10 and Pr = 5.5. For clarity, a coarse-grained vector map 
of size 32 x 32 is shown. The magnitude of the velocity y/ v% + w% is coded in both color and 
the length of the arrows in units of cm/s. The time average is taken over a period of 145 min 
corresponding to 20 000 velocity frames. The two black dots in the figure mark the positions 
(x=-2 1.73cm, z=0cm) (see § 13. 1[) and (x=0cm, z=0cm) (see § 13.2) 1. respectively. 

from the edge. The measured relative temperature differences among eight thermistors 
in the same plate are found to be smaller than 3% of that across the fluid layer. 

During the experiment the entire cell was wrapped by several layers of Styrofoam and 
the cell was tilted by a small angle of about 0.5° (Ahlers, Brown & Nikolaenko 2006) 
so that the measurements were carried out within the vertical plane of the large-scale 
circulation. The mean temperature of water was kept at 29° , corresponding to a Prandtl 
number Pr = v/k = 5.5. The experiment covered the range 5.9 x 10 9 < Ra < 1.1 x 10 11 
of the Rayleigh number Ra = agATH 3 /vn, with g being the gravitational acceleration, 
AT the temperature difference across the fluid layer, and a, v and k being, respectively, 
the thermal expansion coefficient, the kinematic viscosity, and the thermal diffusivity of 
water. As all measurements give the same qualitative results, only that for Ra — 9.5 x 10 10 
will be presented in this paper. At this Ra, the period of the large-scale circulation is 
around 62 sec. 

The details of the particle image velocimetry (PIV) measurements in turbulent RB 
convection have been described and discussed by Xia, Sun & Zhou (2003) and Sun, Xia 
& Tong (2005), here we give only its main features. The laser lightsheet thickness is ~ 2 
mm and the seed particles are 50-/xm-diameter polyamid spheres (density 1.03 g/cm 3 ). 
As the seed particles are neutrally buoyant, they are assumed to follow the motion of the 
fluid. The measuring region has an area of 49 x 49 cm 2 with a spatial resolution of 7.76 
mm, corresponding to 63 x 63 velocity vectors. We chose the vertical rotation plane of 
the large-scale circulation as the laser-illuminated plane, defined as the (x, z) plane. The 
Cartesian coordinate is defined such that the origin (0,0) coincides with the cell center, 
the x axis points to the right, and the z axis points upwards. The horizontal velocity 
component u{x, z) and the vertical one w{x, z) were obtained. The experiment lasted 
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Figure 2. Color-coded contour maps of \uo/u rms \ (a) and \wq / w rma \ (b) obtained at 
Ra — 9.5 x 10 10 and Pr — 5.5. Here, uq and u rms are the time-averaged and the r.m.s. horizontal 
velocities, respectively, and Wo and w rms are the time-averaged and the r.m.s. vertical velocities, 
respectively. 



145 minutes in which a total of 20000 two-dimensional vector maps were acquired with a 
sampling rate of ~ 2.3 Hz. Figure [1] shows the measured mean flow field, together with the 
Cartesian coordinates used in the experiment. In the figure, the magnitude of the mean 
velocity \/u^ + Wq are coded both by color and by the length of the arrows, where uq = 
(u(t)} and wq — (w(t)) are the time-averaged horizontal and vertical components of the 
velocity vector, respectively. One sees clearly that the mean flow is a clockwise rotatory 
motion with a relatively quiet central region and high velocity regions concentrated along 
the perimeter of the cell. To see whether the Taylor's frozen-flow hypothesis is valid or 
not in turbulent RB system, we plot in figure [2] the spatial distributions of the ratios 
between the time-averaged and the r.m.s. velocities, \uo/u rms \ [figure^a)] and \wo/w rms \ 
[figure m 6)], where u rms = yj {[u(t) — uq) 2 } and w rms = yj ([w(t) — wq] 2 ) are the r.m.s. 
horizontal and vertical velocities. If Taylor hypothesis holds, the values of \uo/u rms \ or 
\wo/w rms \ should be much larger than 1. In the figures, it is seen that both \uo/u rms \ 
and \wo/w r ms\ are nearly zero in the cell's central region and \uo/u rm s\ and \wo/w rms \ 
increase, respectively, in the vertical and horizontal directions. The maximization of 
\uo/u r ms\ occurs near the top and bottom plates with a maximum value \uo/u rm s\max — 
2.1, and \wq / w rms \ reaches its maximum value near the cell sidewall with \ wo/w rms \max — 
2.4. For such large r.m.s. velocities, Taylor's frozen-flow hypothesis is then not expected 
to be valid. 

The PIV technique provides us a convenient tool to directly and simultaneously mea- 
sure the local fluctuating velocities at multi-points in a particular plane of interest. With 
the measured u(x,z,t) and w(x, z,t), one can obtain the longitudinal space-time cross- 
correlation functions for both the horizontal and vertical velocities, respectively, defined 
as 

_ {[u(x + r,z,t + r) - u a (x + r, z)][u(x, z,t) -u (x,z)]) 
u u {r,T,z) — — — (z.ij 

UrmsyX ~r *) Z )U rnls yX , ZJ 



and 



([w{x,z + r,t + T) -w (x,z + r)][w(x,z,t) -w (x,z)}) 

C w {r,T\X) = — r r . (2.2) 

Wrrns Z T" » ' J^rms \-E-> % ) 
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Figure 3. (a) Space-time correlations C w (r, r; x) as a function of r and r measured near the cell 
sidewall (a; = —21.73 cm and z = cm). Here, the amplitude of C w (r, t;x) is coded by color. 
(b) The corresponding isocorrelation contours of C w (r,r;x) with the color-coded correlation 
amplitude varying from 0.4 to 0.85 at increments of 0.05 (outer to inner contours). 



3. Results and discussions 

3.1. Near the cell sidewall 

We first study the properties of longitudinal space-time correlations for the vertical ve- 
locity, C w (r, t; x), near the cell sidewall (x — —21.73 cm) at the middle height of the cell 
(z = cm), where the mean and the r.m.s velocities are of the same order [see figure 
QJ&)]. Figure EH a) shows the flood contours of the measured C w (r, r; x) as a function of 
separations r and r and the corresponding isocorrelation contours are plotted in figure 
[3^6). By definition, the maximization of C w {r, r; x) occurs at the origin with a maximum 
value of C w (0, 0; x) = 1. As C w (r, r; x) decays fast near the origin, our present resolution 
could not resolve properly the isocorrelation contours when C w (r,r;x) > 0.85. If the 
Taylor's hypothesis relation (|1.2[) is valid for the present flow field, the isocorrelation 
contours of C w (r, r; x) should be straight lines. However, for the r- and T-range studied, 
one sees in figure [3fb) that the isocorrelation contours of the measured C w (r,r;x) are 
the elongated and closed curves, rather than straight lines, and seem to have a shape 
of elliptical curves that can be described well by (|1.3[) . Furthermore, all isocorrelation 
contours appear to be self-similar, i.e., they share the same preferred orientation and the 
same aspect ratio. One also sees that C w {r, r; x) decays relatively slowly in the preference 
direction, but drops much faster in the direction that is perpendicular to the preference 
direction. 

To determine the characteristic velocities U and V in (|1.3[) . note that from the condi- 
tions drE/dr\ T = and drE/dr\ r = we have, respectively, 

r p = Ut and r p = [U/(U 2 + V 2 )]r, (3.1) 

where r p maximizes C w (r, t;x) for a given r and t p is the pe ak position at which 
C w (r, t; x) reaches its maximum value for a fixed separation r ( He et al. 2OT0j ). Oi 



results for the measured r p as a function of time separation r are shown in figure E5 a). 
One sees clearly that r p increases with increasing r because a longer time lag is needed 
to move velocity fluctuations across a larger separation. It is further seen that the in- 
creasing manner may indeed be described by a simple linear function, r p = Ut, with 
U = 1.52 cm/s. How the measured r p varies with separation r is shown in figure 
Again, one sees that r p increases linearly with increasing r. A linear fit to the data yields 
t p = [U/{U 2 + V 2 )]r with U/(U 2 + V 2 ) = 0.33 s/cm. Note that the similar relation has 
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Figure 4. (a) The measured peak position r p as a function of r. The solid line shows the fitted 
linear function, r p — Ut, with U = 1.52 cm/s. (b) The measured peak position t p as a function of 
r. The solid line shows the fitted linear function, r p = [£//([/ 2 + V 2 )]r, with (7/([/ 2 + V 2 ) = 0.33 
s/cm. All data were obtained near the cell sidewall (a; = —21.73 cm and z — cm). 



also b een observed for the scalar fields in the same system ( Zhou fc Xia 20081 : He et al 



2010h . Taken together, we have U — 1.52 cm/s and V = 1.51 cm/s at x = —21.73 cm 



and z = cm. Comparisons between the fitted velocities U and V and the mean and the 
r.m.s. velocities of the flow will be presented in § 13.31 

With the obtained U and V, we can now test the relation between C w (r,r;x) and 
C w (rE,0;x). Figure [5[a) shows the evolution of C w (r,r;x) as a function of time sep- 
aration r for several different values of r. One sees that the measured C w (r,T;x) all 
have a single peak at the position r p . The peak position r p increases with increasing r 
[see also figure meanwhile the correlation amplitude C w (r,T p ; x) decreases. This 

is because velocity fluctuations at two points decorrelate gradually when the separation 
between these two points increases. For comparison, we first plot in figure the mea- 
sured C w (r,r; x) as a function of the Taylor's separation with calculated from 
(|1.2j) with the mean velocity Uq = 1.95 cm/s, and then shows in figure [5J&) the mea- 
sured C w (r, t; x) as a function of the separation with te calculated from (II .3[) . In the 
figure, only positive parts of C w (r,r;x) are plotted because C w (rE, 0; x) is a symmetric 
function with respect to = 0. It is seen that when using the Taylor's hypothesis [see 
figure [5][ b)] the correlations could not collapse on top of each other, whereas, when using 
the elliptic model [see figure [5[c)] the correlations all collapse well on top of each other, 
indicating that the space-time correlations C w (r,r;x) can be only determined by the 
space correlations C w (rE,0;x) and the solution of ()1.3j) . 



3.2. At the cell center 

Let's now turn to the velocity field at the cell center (x = cm and z = cm), where 
the mean horizontal and vertical velocities are both nearly zero. As the velocity field in 
the cell's central region is approximately locally homogeneous and isotropic (Sun, Zhou 
& Xia 2006; Zhou, Sun & Xia 2008) and we also find that space-time correlations of the 
horizontal and vertical velocities share the same qualitative properties at the cell center, 
only the results of C w (r, r; x) will be presented in this subsection. Figure [6^ a) shows the 
flood contours of space-time correlation C w (r,r; x) as a function of space separation r 
and time separation r and figure ®b) shows the corresponding isocorrelation contours. 
Three features are worthy of note, (i) The measured C w (r, r; x) is a single-peak function 
with the peak locating at the origin and decays with increasing separations r or r. (ii) 
All isocorrelation contours of C w (r, r; x) are closed curves and have an elliptic shape, (iii) 
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Figure 5. The space-time correlations C w (r,r; x) obtained near the cell sidewall (x = —21.73 
cm and z — cm) for various values of r =0, 0.78, 2.33, 3.88, 5.43, 6.98, and 9.31 cm as functions 
of (a) time separation r, (b) the Taylor's separation tt ~ t — Uqt with the mean flow velocity 
Uo = 1-95 cm/s, and (c) the separation tb = \JJr — Ur) 2 + V 2 t 2 with U = 1.52 cm/s and 
V = 1.51 cm/s. 




r(sec.) 



Figure 6. (a) Space-time correlations C' w (r, r; a;) as a function of r and r measured at the cell 
center (x = cm and 2 = cm). Here, the amplitude of C w {r,T\x) is coded by color. (6) The 
corresponding isocorrelation contours of C z (x,r,r) with the color-coded correlation amplitude 
varying from 0.4 to 0.85 at increments of 0.05 (outer to inner contours). 
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Figure 7. (a) The measured peak position r p as a function of r. The solid line marks r p = 0. 
(b) The measured peak position t p as a function of r. The solid line marks r p — 0. All data were 
obtained at the cell center (x — cm and z — cm). 
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Figure 8. (a) The blue solid curves are the isocorrelation contours of Cz{x,r,r) measured at 
the cell center [x = cm and z = cm) with the correlation amplitude varying from 0.4 to 0.85 
at increments of 0.05 (outer to inner contours). Here, we use the same date sets as figure [6] The 
red dashed curves are the elliptic fittings of (|3.2|l to the isocorrelation contours, (b) The fitted 
parameters of the standard elliptic equation (|3.2[) : b vs a. The solid line shows the best fit of 
b = Va, with V = 1.39 cm/s. 



The elliptic isocorrelation contours can be well described by a standard elliptic equation 

2 ? 

t r 

i.e., the isocontours are aligned with the coordinate axis. Here, a and b are two parame- 
ters, determining the lengths of the major and minor axes of the standard ellipse. 

Figure [7^ a) shows the measured r p as a function of time separation r and figure [7^6) 
shows the obtained t p as a function of space separation r. It is seen that r p and r p both 
vary around zero, implying U ~ and [U/(U 2 + V 2 )} ~ [see p.l[) ]. To yield V, we note 
that when [7 = 0, equation (|1.3j) can be rewritten as 

,22 



V 2 t 2 . (3.3) 



Comparing with the standard elliptic equation (|3.2|) . we have b = Va. Therefore, the 
value of V can be estimated from the values of a and b. In figure E^a), we plot again 
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Figure 9. The space-time correlations C m (r,r;x) obtained at the cell center (x — cm and 
z — cm) for various values of r =0, 0.78, 2.33, 3.88, 5.43, 6.98, and 9.31 cm as functions of (a) 
time separation r and (b) the separation t\e = Vr 2 + V 2 t 2 with V = 1.39 cm/s. 



the isocorrelation contours of C w (r,r;x) as the blue solid curves. In the figure, the 
best fittings of (I3.2p to the contours are plotted as the red dashed curves. It is seen 
that the fitted curves collapse well on top of the contours, further confirming that the 
isocorrelation contours have a shape of standard elliptic curves. Figure |8j b) shows the 
lengths of the r-axes of the fitted elliptic curves, b, as a function of the lengths of their 
r-axes, a. The solid line in the figure shows the best fit of b = Va to the data, which 
gives V = 1.39 cm/s. Taken together, we have U ~ cm/s and V = 1.39 cm/s at the 
cell center. 

Figure |9j a) shows the space-time correlations C w (r,r; x) as a function of time sepa- 
ration r for several different values of r. It is seen that C w (r,r;x) obtained at the cell 
center is also a single-peak function. However, unlike the case near the cell sidewall [see 
figure EJa)], the measured peak positions t p here do not vary with r, but all locate at 
positions around r = [see also figure E{£>)]. This is because of the zero mean velocity at 
the cell center. Figure [9)^6) shows the measured C w {r, r; x) as a function of te- One sees 
that when the solution te of (|1.3|) is used, reasonable collapses among these correlations 
are achieved, further confirming the validity of the elliptic model for the present flow. 

3.3. Along the cell's horizontal and vertical central lines 

Figures \TU[a)-(d) show the evolution of the isocorrletion contours of the longitudinal 
space-time correlations C w (r, r; x) measured at the middle height of the cell (z = cm) 
at four different values of x from near the cell sidewall to at the cell center. Figures 
[TCTa) and (d) are the same as figure [3^6) and figure E{6), respectively. We replot these 
figures here for comparison. One sees that the contours are all elliptic closed curves 
with their preferred orientations. Furthermore, the slopes of the preferred orientations 
become smaller as the reference position movies from the wall towards the center. This 
is because the pref erred orientations of the contours ar e directly related to the mean 
velocity of the flow (|He fc Zhandl2006t IZhao fc He! 120091 ) and the mean vertical velocity 
wo, after reaching its m aximum value near the cell sidewall, decreases with the increasing 
distance from the wall ( Qiu fc Tongl200lT ) (see also figured]). 

Direct comparison is made in figure [TT] between the magnitudes of the characteristic 
velocities U (circles) and V (triangles) for C w (r,r;x) along the cc-axis. It is seen that 
\U\ ~ at the cell center and increases along the cell's diameter at the middle height of 
the cell from the cell center to the sidewall and the maximization of \U\ occurs near the 
cell sidewall, while the variation of V is much weaker. Note that the validity of Taylor's 
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Figure 10. The isocorrelation contours of C w (r,r;x) as functions of time separation r and 
space separation r measured at the middle height of the cell (z — cm) at x = —21.73 (a), 
— 13.97 (6), —6.98 (c), and (d) cm. The amplitude of the contours is coded by color and varies 
from 0.4 to 0.85 at increments of 0.05 (outer to inner contours). Note that (a) is the same as 
n gure[3j6) and (d) is the same as figure[6j&). We replot them here for comparison. 



frozen- flow hypothesis requires \U\ 3> V. However, the figure shows clearly that \U\ S> V 
is not the case, i.e., \U\ and V are approximately the same near the sidewall, while at 
the cell center the value of V is even much larger than that of \U\. This further conforms 
that Taylor hypothesis does not hold in the present system. 

To further test the elliptic model, we compare in figure [12] the measured values of the 
characteristic velocities U and V and the ir theoretical predictions Ut and Vt. Based on 
Navier-Stokes equation, Zhao fc Hel ( 20091 ) showed that Ut is a characteristic convection 
velocity proportional to the mean velocity of the flow and Vt is the sum of the random 
sweeping velocity and the shear-induced velocity, i.e. Vt — y / S 2 X 2 + w 2 ms , where the 
subscript "t" indicates theoretical predications and S and A are, respectively, the shear 
rate and the Taylor microscale of the flow. For the vertical velocity profile along the x-axis 
( Qiu fc Tondl200lUSun et ali\2Q05b\ ). the shear rate was evaluated as S ~ 2(Wo) ma x/D 
( He et aO lioiO). where (Wo) max is the maximal value of the vertical velocity along the 
x-axis. The Taylor microscale X(x) was estimated using the equation C w (rE,0;x) ~ 
1 — (r/A(x)) 2 for |r| < 1 cm, where the space autocorrelation function C w (te, 0; x) was 
obtained from the time autocorrelation function C W (Q, r; x) using the relation (11.3[) with 
r = 0. We note that SX(x) <^ w rms (x) for all measuring positions. Comparison is made 
in figure IT2T a) between U (solid circles) and the mean vertical velocity wq (open circles). 
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Figure 11. Comparison of the magnitudes of the characteristic velocities U (circles) and V 
(triangles) for C w (r,T\x) along the x-axis. 
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Figure 12. (a) Comparison of the characteristic velocity U (solid circles) for C w (r, t;x) and 
the mean velocity of the flow wo (open circles) along the x-axis. (b) U vs i»o- The solid line 
shows the linear fit to the data, U = 0.82uio. (c) Comparison of the characteristic velocity V for 
C w (r, t;x) (solid circles) and the theoretical prediction Vt (open circles) of the elliptic model 
along the x-axis. (d) V/Vt vs x. The solid line marks the mean value 1.87 of the ratios. 
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r(sec) r(sec.) 

Figure 13. The isocorrelation contours of C u (r, r; z) as functions of time separation r and space 
separation r measured along the cell's central vertical axis (a; = cm) at z = —21.73 (a), —13.97 
(b), —6.98 (c), and 0(d) cm. The amplitude of the contours is coded by color and varies from 
0.4 to 0.85 at increments of 0.05 (outer to inner contours). 

One sees that both U and wq decrease with increasing x, but the magnitude of U seems 
to be systematically smaller than that of Wq. Figure 1121 6) shows the measured U as a 
function of wo ■ It is seen that U increases with increasing wq and the increasing manner 
may be described by a simple linear function U — 0.82u;o. The solid line in figure IT2T &) 
shows the fitting linear function. Figure [T2Tc) shows the comparison of V and V t . One 
sees that the experimentally measured V is larger than its theoretical predictions Vt for 
all measuring positions. The ratio of V to Vt along the x-axis is plotted in figure [TSTcO. In 
the figure, the dashed line marks the mean value 1.87 of the ratios. One sees that all data 
points vary around the dashed line. This seems to suggest a constant ratio between V 
and Vt, i.e., the experimentally measured V is proportional to its theoretical predications 
V t . 

Finally, we study the longitudinal space-time correlations C u (r, r; z) for the horizontal 
velocity along the cell's central vertical axis. Figures \Wl a)-(d) show the evolution of 
the isocorrelation contours of C u (r, t;z) obtained at x = and four different values 
of z. Again, one sees that all contours are elliptic closed curves with their preferred 
orientations. The preferred orientations guide along the second and fourth quadrants, 
which is due to the negative mean velocities at the measuring positions (see figure [1]), 
and the slopes of the preferred orientations decrease with the increasing distance from 
the plate, which is due to the decrease of the magnitude of the mean horizontal velocity 
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Figure 14. Comparison of the magnitudes of the characteristic velocities U (circles) and V 
(triangles) for C u (r,T;z) along the z-axis. 
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Figure 15. (a) Comparison of the characteristic velocity U (solid circles) for C u (r, t;z) and 
the mean velocity of the flow uq (open circles) along the z-axis. (6) U vs «o- The solid line 
shows the linear fit to the data, U = 0.97u . (c) Comparison of the characteristic velocity V 
for C u (r, r; z) (solid circles) and the theoretical prediction Vt (open circles) of the elliptic model 
along the z-axis. (d) V/V t vs z. The solid line marks the mean value 4.21 of the ratios. 
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(see figure!]}. Figure IT4l shows the comparison of \U\ and V. It is also seen that \U\ and 
V are nearly the same near the top and bottom plates, while V is much larger than \U\ 
at the cell center. The fact that \U\ 3> V does not hold for the horizontal velocity again 
indicates the invalidity of Taylor's frozen-flow hypothesis in the present flow. 

Figure flBT a) shows the comparison of the measured U and the mean horizontal velocity 
Uq along the cell's central vertical axis. It is seen that U and Uq are approximately the 
same for all values of z. Figure HBT 6) shows U as a function of u . The best linear fit to the 
data yields U = 0.97uo ; again indicating that U is proportional to wq- Direct comparison 
is made in figure H^Tc) between V and V t (= \/S 2 \ 2 + u 2 ms ). Here, the shear rate was 
estimated as S ~ 2{Uo) max / H with (Uo) max being the maximal horizontal velocity along 
the z-axis and the Taylor microscale X(z) was evaluated from the space autocorrelation 
function C u (te,0', z). Again, we find that SX(z) <C u rm s(z) for all measuring positions. 
In figure flST c). one sees that similar to the case of C w (r,T;x), the values of V are also 
larger than those of V*. Nevertheless, the ratio V/Vt shown in figure [T5TcO varies around 
its mean value 4.21, suggesting that V is proportional to Vt. Taken together, our results 
reveal that C u (r, r; z) shares the same qualitative properties as C w {r, r; x). 



4. Conclusion 

To conclude, we have presented an systematic experimental study of the velocity field 
in a cylindrical turbulent Rayleigh-Benard (RB) convection cell with unit aspect ratio 
using water as working fluid. The two-dimensional velocity field in the vertical circulation 
plane of the large-scale circulation was measured via the particle image velocimctry 
(PIV) technique and the longitudinal space-time cross-correlation functions for both the 
horizontal and vertical velocities, C u (r, t;z) and C w (t,t;x), were investigated in great 
detail. Our results show that the isocorrelation contours of space-time correlations are 
elliptic closed curves and the space-time correlations C u (r,r;z) and C w (r, t;x) can be 
related to the space correlations C u (rE, 0; z) and C w (rE,0; x), respectively, via the elliptic 
relation (|1.3|) . i.e. r E = (r — Ut) 2 + V 2 t 2 . The characteristic velocities U and V were 
then calculated and studied. We find that the magnitude of U reaches its maximum 
value near the sidewall and plates and decreases when away from the walls and plates, 
while the position-dependence of V is much weaker. Specifically, at the cell center we 
have U ~ and hence the relation (|1.3|) becomes r% = r 2 + V 2 t 2 with its major and 
minor axes coinciding with the r- and r-axis. Note that such relation has the same form 
as Kraichnan's sweeping- velocity hypothesis (|Kraichnanll 1 9641 ) . Direct comparison of the 
values of U and V and their theoretical predictions further show that U is proportional 
to the mean velocity of the flow, while V is systematically larger than its prediction. 

Our results validate the elliptic model in turbulent RB convection, where Taylor's 
frozen-flow hypothesis do es not hold due to the relatively large values of the r.m.s. ve- 
locity. As pointed out by He fc Zhangl ( 2006h . the elliptic model is developed based on 
a second-order approximation, while Taylor's hypothesis implies a first-order approxi- 
mation. Thus, the elliptic model is a generation of Taylor's frozen-flow hypothesis and 
Taylor's hypothesis is only a special case of the elliptic model [i.e., when V — the 
elliptic model's relation (|1.3|) is degenerated to the Taylor's relation (|1.2j) ]. Like Taylor's 
hypothesis, the elliptic model could also be used to translate time series to space series. 
This is because the correlation function is a basic quantity and most statistical proper- 
ties interested in the field of turbulence, such as structure function and power spectrum, 
can be o b taine d theoretically from the correlation functions. In fact, previous work by 
He et all (|201Dh has used the model to translate the temperature power spectrum from 



time domain to space domain. 
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An important implication of the elliptic model is that when r = 0, the relation (|1.3p 
becomes 

r E = (U 2 + V 2 ) 1/2 r. (4.1) 

In this case, r is still proportional to r but just that the proportionality constant is 
(U 2 + V 2 ) 1 / 2 , rather than the mean velocity Uq as stated in the Taylor's hypothesis [see 
(|1.2p with r = 0). This implies that if one is ONLY interested in the scaling exponents 
of C(r, 0) in space domain, such scaling properties can still be correctly obtained by 
studying the scaling of C(0,t) in time domain in the elliptic model even though Tay- 
lor's hypothesis is not valid. We note that Taylor's hypothesis has been widely used to 
study the scaling behaviors of structure functions and power spectru m of the velocity 
and temperature fields in turbulent RB convection ( Lohse fc Xiall2010h . However, it has 



long been known in the field that the condition for Taylor's hypothesis is not often net 
( Lohse fc Xiall201Ch . and hence the results based on Taylor's hypothesis is questionable 



and not convincible. Here, the relation (|4.1j) implies that one does not really need the 
validity of Taylor's hypothesis to reconstruct the space series from the measured time 
series. However, Taylor's hypothesis could not yield the correct scaling ranges in space 
domain, which could only be obtained upon the transform of (|4.1|) in the elliptic model. 



It should be noted that in the original elliptic model (|He fe Zhana 120061 ) the correla- 



tions i?(r, t) = (v(x + r, t + r)v(x,i)) were considered, while we studied the normalized 
quantities C(r, r) = R{r,T)/v rms {x)v rms {x + r) here. The elliptic relation ()1.3[) holds 
for both R(r, r) and C(r, r). However, if the r.m.s. velocity v rms (x) depends on position, 
the characteristic velocities U and V for R(r, r) may be different from those for C(r, r). 
We note that previous works usually focused on C(r, r) (see, e.g.. lHe et al. 2010h . In the 



present work, we have also checked the properties of i?(r, r) and the similar results were 
obtained. 
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